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Abstract 

High-T c superconductivity is unconventional because the gap is not isotropic as 
in simple metals but has d x 2_ y 2 symmetry with lines of nodes. In a fascinating 
thermal transport experiment on a high-T c superconductor, Krishana et al [jl] have 
reported mysterious magnetic field induced first order transitions from a supercon- 
ducting state with gap nodes to a state without gap nodes. We show here that 
this is an experimental manifestation of a novel macroscopic quantum phenomenon 
induced by the magnetic field, qualitatively different from the usual quantum Hall 
effects. It corresponds to the quantization of the superfluid density in a supercon- 
ductor with gap nodes due to the generation of Confined Field Induced Density 
Waves (CFIDW) in the node regions of the Fermi surface. The Landau numbers 
L are not sufficient to index these macroscopic quantum states and the addition 
of a new quantum number £ is necessary. Distinct qualitative implications of this 
non-integer \L, £ > quantization are also evident in a number of recent unexplained 
experimental reports in the cuprates. 
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In conventional superconductors the gap is rather isotropic having s-wave symmetry. 
A breakthrough in the study of high-T c oxide superconductors (HTSC) was the demon- 
stration that instead, the gap is very anisotropic of (i-wave symmetry taking a zero value 
in the so called node directions [Q, []]. Many of the peculiar properties of HTSC are now 
simply understood as resulting from this distinct momentum structure of the gap. How- 
ever, understanding the influence of a perpendicular to the Cu — O planes magnetic field 
on this ground state is actually a fascinating challenge for theorists ||. In a remarkable 
thermal transport experiment, Krishana et al. Jl| have reported that, when a magnetic 
field is applied perpendicular to the CuOi planes, the thermal conductivity shows sharp 
first order transitions from a field dependent regime characteristic of gap nodes 0] to a 
field independent regime indicating the elimination of the nodes. 

The results of stimulated a controversial experimental and theoretical debate which 
was in fact never closed. Thermal transport measurements by other groups || 0, [§] 
have confirmed the surprising findings of Krishana et al.Q]. However, they have also 
reported that the phenomenon is not present in all samples || . Very high quality samples 
are apparently necessary ||. In addition, strong hysteretic behavior was reported in 
Ref. || in disagreement with Refs. |TJ, |7[. Our theory provides a simple explanation 
of this controversy and the sample dependence of the occurence of the phenomenon is 
a major argument for the relevance of our picture. On the theoretical side, the phase- 
transition-point of view has been immediately adopted by Laughlin who suggested that 
this is a transition from d x 2_ y 2 superconductivity (SC) to d x 2_ y 2 + id xy SC ]TT| which is 
nodeless. Models based on the vortex dynamics have also been proposed Jll| explaining 
the saturation of the field effect but could not account for the "kink" and the sharp 
character of the transition. 
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Indications of a field induced transition to a nodeless state are also present in other 
experiments. It has been suggested that the field induced spliting of the zero bias con- 
ductance peak in tunnel measurements may be the result of such a transition HT2| . More 



spectacular are recent measurements by Sonier et al. of the penetration depth in the 
presence of a stronger than usually magnetic field. They also confirm the elimination of 
the nodes by the magnetic field adding a fundamental new element: The elimination of 
the nodes is accompanied by a substantial reduction of the superfluid density [|T3|1 . Such a 
reduction cannot be explained by a phase transition to a new SC state like d x 2_ y 2 + id xy . 

Moreover, a number of seemingly unrelated experimental puzzles emerged recently. 
Neutron scattering experiments reported the generation of AFM moments in the SC 
state by a magnetic field applied perpendicular to the planes [FTI]. The moments appear 



suddenly above a critical field and below a critical temperature [15) . NMR measurements 
have confirmed the phenomenon [jnj. Scanning tunneling microscopy results report in the 
presence of a perpendicular magnetic field a checkerboard structure that covers a region 
around each vortex |T7]]. Finally, measuring both heat capacity and magnetization on an 



extra clean sample of YBa2Cu307, Bouquet et al. |T8| not only confirm the presence of a 
first order transition from a vortex lattice to a so called vortex glass state, but at higher 
fields they observe a surprising transition from the vortex glass state to a new vortex glass 
state which has not found any theoretical explanation so far. 

In this Letter we explore an original physical picture which is shown to account for all 
the experiments. We reveal a new magnetic field induced phase transition from a d x 2_ y 2 
SC state to a state in which d x 2_ y 2 SC coexists with Confined Field Induced (Spin and 
Charge) Density Waves (CFIDW) which develop in the gap node regions of the Fermi 



surface (FS). Field induced density waves have been suggested |L9], [2(], |22|, |23j in 
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order to explain second order metal-insulator transitions in (TMTSF) 2 X (X = PFq, 
CIO4) quasi one-dimensional synthetic compounds under pressure f24]| . The transitions 
we discuss here are novel and the qualitative physics of our CFIDW states is novel as 
well. We consider an hybrid approach. Our HTSC system is built of two subsystems: 
subsystem I is the Fermi Surface (FS) region covered by the superconducting gap and 
subsystem II is a virtual normal quasiparticle region created by the magnetic field and 
centered in the node points of the FS. The CFIDW will eventually develop in region II 
because of the orbital effect of the field in this region. The orbital effect of the field in 
region I induces vortices that are assumed to be irrelevant. The Zeeman effect is not 
considered here for simplicity, however in region II it can be shown to further stabilize 
our CFIDW states while in region I it is negligible. The relative momentum extension of 
regions I and II is unknown and fixed by the energetic competition of CFIDW with SC. 
The above picture has some similarity with the partially depaired Fulde-Ferrell-Larkin- 



Ovchinikov state |2^, where the normal state is created by the magnetic field over a 
portion of the Fermi surface and competes with superconductivity. 

In HTSC the gap is d x 2_ y 2 with nodes in the (±7r, ±7r) directions where region II is 
centered. In region II we necessarily have open FS sheets. Therefore, we can write the 
dispersion of subsystem II in the form: ££ J = hp) — 2t 2 cos(k 2 /X) — 2t 2 cos(2k 2 /X) 

where X is the unknown momentum extension of region II (see Fig. la). k\ is along 
the (±7r, ±7r) directions perpendicular to the open FS sheets of subsystem II, and k 2 
is perpendicular to k\ and therefore along the open FS sheets where we keep only two 
harmonics without influence on the generality of the results. 

We are able to calculate explicitely the spin and charge susceptibility of subsystem 
II in the presence of the magnetic field exploiting methods developed for the study of 
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(TMTSF) 2 X compounds p0| , 21| . The details of the calculation and a discussion with 



various electronic dispersions in II will be given elsewhere. We report here the simplest 
results sufficient to illustrate the surprising physics of the system. Constraining the nesting 
in the (±7r, ±7r) directions one can show that a first order field induced density wave gap 
in region II is given by A DW = W exp{ — [gN(Ep)I'l(X)]~ 1 } where 

U 2 X \ ( 2t' 2 X 



)pj \eHvp, 

Here J n {x) are Bessel functions, L is the index of the Landau level configuration, e is the 
charge of the electron, H the magnetic field, N(Ep) the density of states at the Fermi 
level (in region II), g a scattering amplitude of Coulombic or phononic origin, W the 
bandwidth in the (jr, it) direction and Vp the Fermi velocity. Higher order gaps are not 
reported here for sake of clarity. 

CFIDW states will develop only if the absolute free energy gain of the system due to 
the opening of a CFIDW gap in region II is bigger than that lost by the elimination of 
SC from this region. This condition implies the following inequality: 



l£ Jfe F sin(#7r/2) > 
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In 



2W 2 nZ 



A2 c [7rZ-sin(7rZ)]J 



-i 



(2) 



gN(E F ) 

Z is a dimensionless number varying from to 1 that represents the relative extension 
of the CFIDW states over the FS and is related to X by X « kpsm(Zir/2). A sc is the 
maximum value of the d x 2_ y 2 SC gap. The CFIDW state must be confined in momentum 
space with a DW gap smaller or equal to the absolute superconducting gap in the borders 
of region II. We therefore have: 

W 



I 2 L (k F sm(Zn/2)J < 



gN(E F )\n- 



(3) 



A sc sin(iv7r/2 

The equality in (3) fixes the relative extension Z of the CFIDW (for each L) and 



therefore I\ ykp sin(Z7r/2)J which then fixes the CFIDW gap Adw and the critical tem- 
perature Tdw at which the CFIDW forms. A graphic solution for Z is shown in Fig. lb. 
Surprisingly, there are two possible values of Z for a given Landau level configuration. 
Therefore, the Landau numbers L are not sufficient to index our quantum states. A new 
quantum number (, associated with the two possible momentum extensions Z for each L 
configuration must be added. Physically ( will index the quantization of the superfluid 
density. In fact, each quantized value of Z corresponds to a different relative extension of 
the SC region over the FS and therefore to a different density of superfluid carriers. 

Using for our parameters values extracted from the experiments on HTSC and assum- 
ing a conventional scattering gN(Ep) « 1 we obtain results like those reported in Figures 
2 and 3 in remarkable agreement with the data of []TJ and |Tj|. In Fig. 2 is reported 
the dependence of T DW on the magnetic field in the various \L,( > configurations. We 
also plot in this figure the experimental points of refs [Q] and [I3|. In Fig. 3 we plot 
the corresponding magnetic field dependence of the accessible relative extensions Z of the 
CFIDW in each \L, £ > configuration for the same parameters. A quantitative fit of the 
experiments as in figure 2 establishes that the orders of magnitude of the involved param- 
eters are compatible with the experimental data. More importantly, there are distinct 
qualitative experimental facts which strongly support our picture. As one can see in Fig. 
2, the Tdw versus critical magnetic field profile of the data in |I| (circles) show a "repta- 
tion" shape. The first two points have a bigger field slope than the next two points and 
so on. To the best of our knowledge, no explanation of the reptation behavior has been 
reported so far. Within our analysis this "reptation" profile is due to the quantization of 
the CFIDW states. Each slope corresponds to a different \L,( > quantum configuration 
of the system. The higher the field is, the smaller is the field slope in agreement with the 
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experiment. Moreover, if one associates the vortex solid to vortex glass transition with the 
SC to SC+CFIDW transition, one naturally understands the unexplained vortex glass I to 



vortex glass II transition of Bouquet et al jlfl/ as a transition between two different \L, £ > 
configurations of the SC+CFIDW state. The higher field \L,( > state corresponds to a 
higher relative extension of the CFID W state and therefore a less rigid vortex structure 
as is experimentally observed. Alltogether, this provides a microscopic viewpoint for the 
origin of the spin-glass states. 

The remarkable sensitivity of the occurence of this phenomenon on sample quality 
||, further corroborates our picture. Within our analysis, a small magnetic field affects 
the large SC gap of HTSC creating CFIDW's in the node region because of momentum 
confinement which constrains the cyclotron orbits to be large in real space and conse- 
quently the involved flux is large as well. The occurrence of the CFIDW's requires very 
clean samples because the mean free path must be bigger than the cyclotron orbits. Sam- 
ples which may appear of high quality from the usual criteria (width of the SC transition 
or magnitude of the SC T c ) may not be sufficiently clean to support large cyclotron orbits 
and show the CFIDW states. This explains the puzzling sample dependence of the oc- 
curence of the phenomenon. The minimum field at which the phenomenon is observable 
is also limited by sample defects because if fields are too small, momentum extensions 
of the CFIDW are also too small (see figure 3) and consequently the required cyclotron 
orbits are too large. Our d x 2_ y 2 to d x i_ y i + CFIDW transition appears only above 1 
Tesla in |5| while it is already present at 0.6 Tesla in |TJ because samples in [I] are cleaner 
admiting bigger cyclotron orbits. Because all \L,C > configurations are nearly degenerate 
for the total system, the occurence of magnetic hysteresis on a dynamic probe like ther- 
mal transport will depend on the exact conditions of the magnetic cycle and on sample 
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quality as well. In the magnetic cycle of Ref. |J showing hysteresis, the field orientation 
is reversed when the field maxima are reached (the field passes through zero), while in 
Refs. |l this is not the case and hysteresis is absent. 

Our analysis is also the first to establish a natural relationship between the thermal 
transport [0 and penetration depth Jl3| puzzles. In |13| CFIDW develop at a much 
higher T^w for a given field than in [I]] because the system is cooled in the presence of 
the magnetic field. By field cooling the sample, the first accessible \L,( > configuration 
is the one with the higher Tow which also corresponds to the higher Z (cf. figure 3). On 
the other hand, in the experiment of Jl|, the temperature is kept constant when the field 
varies. Not only we reproduce simultaneously the experimental T^y/ versus field profiles of 
both |fj] and [I"3|, but we also account for the reduction of the superfluid density reported 



in [pl|. As one can see in figure 3, in the regime 4 - 6 T explored in [13 , the CFIDW 
occupy about a quarter of the FS [Z pa 0.25 to 0.35) which is in very good agreement 
with the ~ 15% to 25% reduction of the superfluid density reported in |13|]. Furthermore, 
our CFIDW states correspond to a real space pattern that has all the characteristics of 



the one observed by STM [l7j. 

Finally, in agreement with all the experiments, our CFIDW states appear only for 
fields perpendicular to the planes because for fields parallel to the planes cyclotron effects 
are irrelevant. 
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Figure Captions 



Figure 1: a): Schematic view of the SC-CFIDW competition, b): Graphic solution 
of (3) in the L = 1 configuration and a field of 3 Tesla. Two different relative extensions 
Z (JT ~ kp sin(Zvr/2)) of the CFIDW are possible. 

Figure 2: (color): Critical temperature Tpy/ versus critical field for the formation of 
a CFIDW state in the different quantum configurations \L, ( >: L = 4 (orange), L = 3 
(red), L = 2 (blue), and L = 1 (green). In all L configurations, full lines and dotted lines 
correspond to the two different ( configurations (full lines to the higher Z solution). The 
L = lines are not shown for clarity. The open circles are the corresponding experimental 
points of PJ and the open squares the experimental points of ||13||. 



Figure 3: (color): Relative extensions Z versus magnetic field for various quantum 
configurations \L,( >■ The different colors and line-styles correspond to the same \L,( > 
as in Fig. 2. 



13 



Magnetic Field (Tesla) 

O K> -P^ ON 




